Spin Dynamics of a J1-J2-K Model for the Paramagnetic Phase of Iron Pnictides 
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We study the finite-temperature spin dynamics of the paramagnetic phase of iron pnictides within an antifer- 
romagnetic Ji — J2 Heisenberg model on a square lattice with a biquadratic coupling —K{Si ■ Sj)"^ between 
the nearest-neighbor spins. Our focus is on the paramagnetic phase in the parameter regime of this Ji — J2 — K 
model where the ground state is a (tt, 0) coUinear antiferromagnet. We treat the biquadratic interaction via 
a Hubbard- Stratonovich decomposition, and study the resulting effective quadratic-coupling model using both 
MSW and SBMF theories; the results for the spin dynamics derived from the two methods are very similar We 
show that the spectral weight of dynamical structure factor 5(q, co) is peaked at ellipses in the momentum space 
at low excitation energies. With increasing energy, the elliptic features expand towards the zone boundary, and 
gradually split into two parts, forming a pattern around (tt, n). Finally, the spectral weight is anisotropic, being 
larger along the major axis of the ellipse than along its minor axis. These characteristics of the dynamical struc- 
ture factor are consistent with the recent measurements of the inelastic neutron scattering spectra on BaFe2As2 
and SrFe2As2. 
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I. INTRODUCTION 

The emergence of superconductivity in iron pnictidesi"^, 
near an antiferromagnetically ordered state^^ in the phase di- 
agram suggests strong interplay between the superconductiv- 
ity and magnetism in these materials. Elucidating the mag- 
netic excitations is therefore important for understanding not 
only the overall microscopic physics of these systems but 
also their superconductivity. In the parent compounds, the 
observed(7r, 0) antiferromagnetic order arises either within a 
weak-coupling approach invoking a Fermi surface nesting.^^ 
or from a strong-coupling approach whose starting point is a 
local moment Ji — J2 — K model^"— 

The strong-coupling approach is based on the proximity of 
the metallic ground state of the parent pnictides to a Mott lo- 
calization transition, which gives rise to quasi-local magnetic 



moments 
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This incipient Mott picture corresponds to a 
ratio of U (a measure of the Coulomb repulsions and Hund's 
couplings among the Fe 3d electrons) to t (the characteristic 
bandwidth of the Fe 3d electrons) which is not too far be- 
low the Mott threshold Uc/t, which is usually of order unity. 
This is supported by many experimental observations. For 
instance, the room-temperature electrical resistivity of par- 
ent iron pnictides is so large (even when the residual re- 
sistivity is relatively small signaling the smallness of elas- 
tic scattering) that the extracted mean-free path of quasipar- 
ticles would be comparable to the Fermi wavelength; this 
is typical of bad metals near a Mott trasition. Similarly, 
the Drude weight in optical conductivityi^ii^ is strongly sup- 
pressed from its non-interacting counterpart, providing a di- 
rect measure of the proximity to the Mott transition. This 
is further corroborated by the the temperature-induced spec- 
tral weight transferrin— which is also characteristic of met- 
als near a Mott transition. In the spin sector, zone boundary 
spin waves have been observed by inelastic neutron scattering 
(INS) measurements in the magnetically ordered state of sev- 
eral 122 iron pnictides compounds.^' Both the large spectral 
weight and the relatively-small spin damping suggest quasi- 
localized moments, which are expected near the Mott transi- 



tion where the spin excitations arise out of incoherent elec- 
tronic excitations. Additional evidence for the incipient Mott 
transition picture has come from the observation of a Mott 
insulating phase in the iron oxychalcogenidesi^ This mate- 
rial contains an expanded Fe square lattice compared to the 
iron pnictides, which reduces t, thereby enhancing U/t be- 
yond Uc/t (Ref. I22I). Likewise, the Mott insulating behavior 
of the alkaline iron selenides^J can also be interpreted as the 
result of a reduced effective t and, correspondingly, an en- 
hanced U/t beyond Ujti^ 

In the vicinity of Uc/t, where correlations are strong, it is 
natural that the spin Hamiltonian contains not only two-spin 
interactions, such as Ji and J2 Heisenberg exchange between 
nearest- and next-nearest- neighbor spins on a square lattice, 
but also interactions involving higher number of spins. These 
naturally include, for instance, the ring-exchange coupling in- 
volving four spins on a plaquette, and the biquadratic coupling 
of the form —K{Si ■ S^)^ in systems with spin size S > 
The subject of the present study is to show how such non- 
Heisenberg interactions, particularly the biquadratic interac- 
tion, influence the spin dynamics in the paramagnetic phase. 

Spin dynamics in the parent iron pnictides have been most 
extensively studied in the low-temperature state (T < T/v) 
with both antiferromagnetic order and orthorhombic structural 
distortion. Here, the INS experiments up to high energies 
(on the order of 200 meV) show that the spin wave excita- 
tions in these compounds are highly anisotropic, with a dis- 
persion which can be understood in terms of an anisotropic 
•hx — Jiy — J2 model with Jix ^ Jiy 1^^*212^. The anisotropy 
in the nearest-neighbor coupling is compatible with the or- 
thorhombic structure, and its degree could reflect an orbital 
orderingisSi^ Detailed theoretical studies of the magnetic ex- 
citations in the ordered phase have been carried out in such a 
•hx — Jiy — J2 model^ and in a Ji — J2 — if modeli^i^ 
It should also be noted that terms such as the biquadratic cou- 
pling could be inferred from the sublattice angle dependence 
of the ground-state energy in LSDA calculation^, and were 
shown to appear naturally as a result of the orbital ordering 



between Fe A^z and d^^ orbital: 
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Our focus is instead on the spin dynamics in the paramag- 
netic phase of the parent iron pnictides, which has only re- 
cently been studied experimentally. The initial work by Di- 
allo et al.^- measured the spin dynamics of CaFe2As2 at rela- 
tively low energies, below 70 meV. Theoretically, four of us^ 
studied the spin dynamics in the paramagnetic phase of the 
Ji — J2 model (with or without an additional fermion damp- 
ing). We showed that the experimentally observed elliptical 
features of the spin spectral weight in momentum space are 
well-described by this model and we determined the change 
to the elliptical features at high energies. 

More recently, Harriger et ai— reported measurements of 
the spin dynamics in the paramagnetic phase up to high en- 
ergies (above 200 meV) in BaFe2As2. The INS measure- 
ments confirmed the quasi-two-dimensional spin dynamics 
found at low energieSj^ and characterized the evolution of the 
low-energy elliptic features as they expand towards the zone 
boundary as the energy is raised, and determined the high- 
energy dispersion which appears to require a Jix ^ Jiy de- 
scription even though the paramagnetic phase has a tetrago- 
nal structure. Similar data have also been reported by Ew- 
ings et al. in SrFe2As2 42 Theoretically, Park et alM- analyzed 
the spin dynamics in the paramagnetic state within a dynam- 
ical mean-field theory (DMFT) for interactions U/t ^ Uc/t, 
demonstrating that the DMFT approach captures key features 
of the neutron scattering results, including the ellipticity of the 
map of the structure-factor peak in the Brillouin zone. 

In this paper, we study the spin dynamics of the Ji — J2 — K 
model in the tetragonal paramagnetic phase using both mod- 
ified spin wave (MSW) and Schwinger boson mean-field 
(SBMF) theories. The results from the two methods are in 
very good quantitative agreement with each other We show 
that, for a moderate biquadratic coupling K, the dynamical 
structure factor S{q,uj) has not only elliptic features near 
(tt, 0), which expand with increasing energy and split into 
peaks surrounding (tt, tt), but also an anisotropic distribution 
of the spectral weight that is larger along the major axis of 
each ellipse than along its minor axis. These properties agree 
well with the INS experiments^^iili 

The remainder of the paper is organized as follows. In 
SeclUwe introduce the Ji — J2 — K model and describe the 
MSW and SBMF theories used in this paper In Sec. Ulllwe 
show how the biquadratic coupling K influences the mean- 
field phase diagram and magnetic excitation spectrum. In 
Sec. |lV]we calculate the dynamical structure factor S{q,uj). 
We also show that the spectral weight exhibits anisotropic fea- 
tures, discuss the evolution of the anisotropic features with in- 
creasing excitation energy, and explain how these properties 
arise from our theory. In Sec.|V]we first discuss some possi- 
ble generalizations of the Ji — J2 — K model we are studying 
in this paper In the same section, we then consider the effect 
of itinerant electrons, and compare our study with other theo- 
retical approaches to the spin dynamics. Sec.|VT]is devoted to 
a comparison with the INS experiments on the paramagnetic 
phases of the parent 122 iron pnictides, and Sec I VIII contains 
a few concluding remarks. In three appendices, we expound 
on the Ising transition at small Ji / J2 ratios, and discuss the 
effects of both the ring-exchange interactions and interlayer 



exchange couplings. 

II. MODEL AND METHODS 

The Ji — J2 — K model is defined on a two-dimensional 
(2D) square lattice with the following Hamiltonian: 

H = Ji ^ Si ■ Si+s + J2 ^ Si • Si^s' 

i.S i.S' 

-KY,{S^■S^+S)\ (1) 
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where Ji and J2 respectively denote the antiferromagnetic ex- 
change couplings between spins located in the nearest neigh- 
bor [5 = x,y) and next-nearest neighbor (6' = x zk y) sites. 
K is the coupling for the biquadratic interaction between the 
nearest neighbor spin pairs. 

To fully explain the experimentally observed (tt, 0, tt) an- 
tiferromagnetic order, an exchange coupling along the third 
dimension, should also be included. However, we find the 
model defined in Eq. ^ already allows us to understand the 
experimentally observed quasi-2D spin dynamics. Hence, we 
concentrate on this 2D model in the main text, and discuss the 
influence of the interlayer coupling Jr on the spin dynamics 
in AppendixICl 

The Hamiltonian of Eq. ([T]i is studied using both MS Wi^^dl 
and SBMF"*^ methods. Here, we focus on the parameter 
regime where the ground state has a collinear (tt, 0) anti- 
ferromagnetic order, and decompose the biquadratic interac- 
tion term of the Hamiltonian using two Hubbard-Stratonovich 
fields The effective Hamiltonian reads as 

H = Ji ^ Si • S,+5 + J2 ^ Si • Si^s' 

i.S i.S' 

-2A'^r,,5S,-s,+5 + x^r2,. (2) 

i.S i.S 

At the mean-field level, the Hubbard-Stratonovich fields are 
treated as static quantities, and can be expressed using equal- 
time spin correlators as: Ti s = (S^ • Si+s)- The Hubbard- 
Stratonovich transformation itself is exact. The static approxi- 
mation is made in accordance with the level of approximation 
inherent to the MSW and SBMF methods, which incorporate 
static self-energies for the respective boson fields. As shown 
below, our approach has two important features: i) it is capa- 
ble of studying the Ising correlations at nonzero temperatures; 
and ii) the MSW and SBMF approaches yield consistent re- 
sults. 



A. The modified spin wave theory 

The MSW theory^2il^ has been applied to the Ji — J2 model 
by four of usM In this approach, a local spin quantization axis 
is defined at each site along the classical ordering direction 
nf. The Hamiltonian in Eq. ^ is then expressed in terms of 
Dyson-Maleev (DM) bosons via a local DM transformation: 
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Si-nf = S~ ajoi, S+ = V2S{1 - a\a^/2S)a^, and S" = 
\/2Sa\. Minimizing the free energy under the constraint of 
zero sublattice magnetization {S — a\ai) = by introducing 
a Lagrange muhipher fjL, and with respect to (= Vi.s, by 
assuming translational symmetry), we obtain 



• 2 r 2 



„2 V^^2 

2' 



(3) 



where 4> = arccos(Jlf • ^f^^). fs = (alai+s) and gs = 
{aiUi+s) are the ferromagnetic and antiferromagnetic bond 
operators, respectively. Minimizing the free energy with re- 
spect to gives either sin = for nonzero fs and gs, or (f> 
can be arbitrary if fs = gs = 0. This defines two phases sep- 
arated by a mean-field temperature scaleifi Tcro^ at T > T^ro, 
(f> is arbitrary, and the system has Cav lattice rotational sym- 
metry; while for T < T^o, the €4^ symmetry is broken and 
the system is Ising ordered, corresponding to either = or 
^ = TT. In MSW theory, the Ising order parameter can be de- 
fined as a = 2(cos2 + g^) - sin^ + gl)). From 
Eq. (O, if we define r± = (Tx ± Ty)/2 as the symmetric 
and antisymmetric Hubbard-Stratonovich fields, we find that 

r_ = cr/4. 

Minimizing the free energy with respect to other variational 
parameters, we obtain a set of self-consistent equations: 



fs = niQ + 



gs = mo + 



k 
k 



i ) cos(k • S), (4) 



cos(k-5'): (5) 



1 



k 



(6) 



where J\f is the total number of lattice sites, S 
S' = x,y,i ± y. In Eqs. (HI-®, 

Ak = 2 sin^ ^Jixgx cos + 2 cos^ '^Jiy9y cos ky 

+'iJ29x+y cos k^ cos ky, (7) 

Bk = 'iJ29x+y - A* + 2sin2 t:iJix9x - JlyfyC^ - cosky)) 



+2cos^ -{Jiy9y - Jixfx{l - cosk^)), 



(8) 



and the Bogoliubov angle 6'k is defined via tanh26'k = 
^k/i^k- The boson dispersion ek = y/ B^ — and the bo- 
son number rik = [Gxp(£k/fcs7') — 1]^^. At T = 0, the spec- 
trum of the DM bosons becomes gapless at wave vector Q and 
0. This corresponds to a long-range antiferromagnetic order 
at Q 7^ with a nonzero spontaneous magnetization mp. In 
this case, the summation J2k runs over all k values that make 
Ek > 0, and the contribution from the = terms is taken 
into account separately by mo. For T > 0, niQ = 0, and the 
system is paramagnetic. Here the summation is performed in 
the full momentum space. In the presence of a small third- 
dimension coupling Jj, there will be a nonzero mean-field 
Neel temperature, Tjvo; this is discussed in AppendixICl 



Note that these self-consistent equations are exactly the 
same as those for the isotropic Ji — J2 model.^"^ But the defi- 
nitions of Ak and i?k are different. In Eqs. (|7]i and (H)) above, 
we defined the effective exchange couplings Ji^ (Jiy) along 
the x{y) direction as follows: 



Ji 



x{y) 



Ji - 2A"r 



x{y), 



(9) 



expressed in terms of the Hubbard-Stratonovich fields r^(j^) 
of spin-spin correlators in Eq. ([3]). Although in the J1 — J2 — K 
model the bare nearest neighbor exchange coupling Ji is 
still isotropic, a nonzero biquadratic coupling K leads to an 
anisotropic effective coupling Jix 7^ Jiy in the Ising ordered 
phase where 7^ Tj^, i.e. the nearest-neighbor spin correla- 
tors along X and y are unequal, similarly to the situation found 
originally^ for the Ji — J2 model. 



B. The Schwinger boson theory 

In the Schwinger boson representation^^ the SU(2) spin op- 
erators are rewritten in terms of two Schwinger bosons via 
the transformation: = ^{alui — blbi), = albi, and 
Sj^ = bjui. To limit the boson Hilbert space to the physi- 
cal sector, a constraint ajoi + blbi = 2S is imposed on each 
site. This can be generalized to the case of either SU(-/V)^ 
or SP{N)^^^ spins, in either case there will be N boson de- 
grees of freedom at each site. For the experimentally observed 
Q = (tt, 0) or (0, tt) antiferromagnetic collinear phase in the 
122 parent compounds, the (a6)-plane spin-spin correlations 
are expressed as: 



+e(*,j)[2/t.4-- 5(5+1)], 



(10) 



X, y, and where fs = fi,i+s = ^(a|a,+6 + bjb^+s) and gs = gis+s = 
i {ttibi-^s ~ biGi-^s) are respectively the ferromagnetic and an- 
tiferromagnetic bond operators. The function = 1 if 
i and j are on the same stripe sublattice, and 6(i, j) = 
if i and j are on different stripe sublattices. The Hubbard- 
Stratonovich field is then Ts = — \gs\'^, and in the case 
of (tt, 0) ordering we find 



-9I, 



— fy- 



(11) 



Comparing this to Eq. Q, we see that the spin correlators 
coincide with those in the MSW theory if one sets </) = tt. 
Similarly, the case of (0, tt) ordering corresponds to = in 
Eq. (O. In both cases, and Ty have opposite sign, leading 
to the anisotropy in the effective spin-spin exchange couplings 
Jix 7^ Jiy, from Eq. 

By introducing Fourier transformation^ 



77^ 



1 

T 



awe 



^6ke 



:(k-2) 



:(k+S)., 



(12) 



(13) 
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and making a Bogoliubov transformation to a new quasi- 
particle creation/annihilation operators ~ cosh6'kak + 
isinh0fc6^j., one arrives at the mean-field free energy den- 
sity, which can be generaUzed to the Sp(A^) formii 



NT ^ r / cjk 



iVA 



8 



.961 



(14) 



where z is the coordination number, and A is the Lagrange 
multiplier associated with the imposed constraint that on av- 



erage, the number of bosons per site X]^=i '^icr = NS. Here 
Wk — (-Bk — A)2 — A? is the dispersion of the Bogoliubov 
quasiparticles, expressecpi in terms of the variables 



i E Js gs e" 



ikS. 



Bu = ^Js Is e" 



(15) 



the BogoUubov angle tanh20k = ^k/(A — i?k)- The dis- 
persion relation Wk explicitly depends on the ordering wave- 
vector Q and has minima around k = ±Q/2. In the regime 
when J2 > Ji/2, the minimization of the free energy results 
in Q = (ti'iO) or (0,7r). For example, for Q = (tt, 0), the 
expressions for ^k and Sk become: 

Ak = 2JixgxSmkx + 4 J25a;+y sin fc^; cos fcy , (16) 

Bk = 2Jly/y COSfcy. (17) 

In the large-iV limit of the Sp(A^) spin, the mean-field free 
energy Eq. ( fT4l i becomes exact4iii The observable magnetic 
excitation spectrum is obtained from Wk by a Q/2 shift: £k = 
'^k-Q/2- At T = 0, the magnetic order results in the gapless 
Goldstone modes at k = and Q, as expected. The SBMF 
theory is known to reproduce well the spectrum of spin waves 
in both ferro- and antiferro-magnets.i^^ 

Below, we focus on the paramagnetic phase at T > 0, with 
short-range Q = (tt, 0) antiferromagnetic correlations (the 
case Q = (0, tt) is obtained by C4 lattice rotation). We obtain 
the following self-consistent equations from the saddle-point 
minimization of the free energy Eq. ( fT4l i: 



k 



A 



"k 



i)cos(k.5), (18) 



1 



gs' = 



^ cjk 

k 



n-k 



sin (k • 5') , 



S 



Sk-A 



Wk 



(19) 



(20) 



where S = y, and 6' = x, x zt y. Under the transformation 
Bk - A ^ Bk and k ^ k - Q/2, Eqs. in the 

SBMF theory and Eqs. ©-([Hll in the MSW mean-field theory 
have exactly the same form in the short-range (tt, 0) correlated 
paramagnetic phase. Therefore, the two methods yield exactly 
the same mean-field phase diagram and boson dispersion, as 
corroborated by explicit numerical comparison. We further 
verified that these two theories give similar results for the spin 
dynamics of the Ji — J2 — K model. 




FIG. 1 . (Color online) Mean-field phase diagram in the MSW theory 
for S — 1 and J1/J2 = 1. Phases I, II, and III respectively denote 
the (tt, 0)/(0, tt) long-range antiferromagnetically ordered state (at 
T = 0), the Ising ordered paramagnetic state, and the isotropic para- 
magnetic state. The solid red curve refers to the mean-field temper- 
ature scale Tcto = To- In the shaded region the effective exchange 
coupling Jiy < 0. 



III. MEAN-FIELD PHASE DIAGRAM AND EXCITATION 
SPECTRUM 



Since INS measurements suggest Ji ^ J2 for several 122 
compoundS)2ii^2i^ our discussion on the Ji — J2 — K model 
is focused on this parameter regime. Fig. [T] shows the mean- 
field phase diagram of the 2D Ji — J2 — K model using the 
MSW method for S* = 1 and Ji/,/2 = 1- We identify three 
different phases. Phase I coiTesponds to the (tt, 0)/(0, tt) an- 
tifeiTomagneticaly long-range ordered phase; it exists only at 
T = in the 2D model. Phase II and phase III are both para- 
magnetic. They are separated by a mean-field Ising transition 
temperature T^q. We find that for J1/J2 = 1, this transi- 
tion is first-order, as shown in Fig. [J] But it can be either 
first-order or second-order for J1/./2 0-9' as discussed in 
more detail in Appendix [A] In the low-temperature phase II, 
either fx 7^ fy or 7^ gy (see Fig. |2]i, corresponding to 
an Ising ordered phase with either (tt, 0) or (0, tt) short-range 
antiferromagnetic correlations. This Ising ordered phase al- 
ready exists in the isotropic Ji — J2 modeli^MLiS But here, we 
find that a nonzero K enhances T^o, and K drives the effec- 
tive nearest-neighbor exchange couplings to be anisotropic. 
As shown in Fig. |2l in the (tt, 0) Ising ordered phase (cor- 
responds to <f> ~ it), the effective coupling Jiy can even be 
ferromagnetic. This is important for understanding the exper- 
imentally observed anisotropic magnetic excitations at high 
energies in Ca-123^ and Ba-122.^^ Phase III at T > T^q is 
the Ising disordered paramagnetic phase. In this phase the 
effective nearest-neighbor exchange couplings are isotropic 
because the nearest-neighbor bond correlators are zero. But 
the next-nearest-neighbor bond correlations may still be fi- 
nite in this phase. One may define another temperature scale 
To, above which the next-nearest-neighbor bond correlations 
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FIG. 2. (Color online) The temperature evolution of the mean-field 
parameters in the MSW theory for S — 1, Ji j Ji — 1, and K/ J2 = 
0.8. 



vanish and the system are decoupled into isolated local mo- 
ments. Note that Tq does not refer to a phase transition, and 
the discontinuity of the bond correlations at Tq is an artifact 
of the mean-field theory 4^ In general, Tq and T^o are two dif- 
ferent temperature scales satisfying Tq ^ T^rQi^^^ But for 
J1/J2 ^ 0.9, Tq = Tcto for any K/J2 ratio, as shown in 
Fig. ID The phase diagram obtained in the SBMF theory is 
identical to the one shown in Fig.[T] 

The finite coupling K not only changes the phase boundary 
of the mean-field phase diagram, but can also dramatically 
influence the boson excitation spectrum. In Fig.[3ja), we show 
the dispersions of the DM and Schwinger bosons along two 
high-symmetry directions in momentum space for various K 
values in phase II with (j) ~ tt using the same parameters as in 
Fig.[T] We see that the dispersion in Schwinger boson theory 
matches the one in the MSW theory exactly. 

The dispersion shows a gap at (tt, 0) (and also at (0, 0)), 
with the size 



x+y 



(21) 



At low temperatures the gap is small since /i as T <C 
T^Q. In this limit, the excitation near (tt, 0) can be approx- 

where the 



imated by Sk = y'^ixi^ 
velocities are respectively: 

'iJ29x+y + '2Jlx9x, 



A?, 



(22) 



"la = V {^-hgx+y + '2.Jlxgx){^J2gx+y - '2-Jlyfy) + 2Jiyfy^. 

(23) 




(0,0) 



(K,0) 



(71,71) 



200 



> 



100 




(0,0) 



(71,0) 



(7t,7l) 



FIG. 3. (Color online) (a): MSW dispersion along high symmetry 
directions in the paramagnetic Brillouin zone for the 2D Ji — J2 — 
model at S = 1, J1/J2 = 1.0, T/J2 = 1.0, and various K values. 
For comparison, the red dashed curve shows the dispersion in the 
SBMF theory for the same parameters and K/J2 = 0.8. The gaps 
at (0, 0) and (tt, 0) are too small to be seen in the figure, (b): The 
symbols show the dispersion from the INS data at T = 150 K in 
BaFe2As2, taken from Ref. The data can be fit by any of the 
theoretical dispersion curves [as in (a)] that lie within the shaded 
region. 



The excitation develops to the gapless Goldstone mode at T = 
0. At (tt, tt) (and also at (0, tt)) the dispersion has a different 
gap 



A. 



{^■hgx+y - ^Jlyfy ^ fJ-){4^Jlxgx " 4Jiy/y 



m). 

(24) 

The features that vix ^ I'ly and Ai 7^ A2 already exists in 
the isotropic Ji — J2 model. In the Ji — J2 — K model, Ai 
at (tt, 0) is only weakly affected by K because it is dominated 
by /.J. But A2 at (tt, tt) is strongly influenced. It increases 
with K. For sufficiently large K, approximately where Jiy 
changes sign to be ferromagnetic (the shaded region in Fig.[T]i, 
the dispersion at (tt, tt) turns from a local minimum to a maxi- 
mum, as shown in Fig.jSja). Similar behavior in the spin-wave 
dispersion of the Ji — J2 — K model has also been discussed 
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form in both MSW and Schwinger boson theories. 



(c) 



-1 



-2 





FIG. 4. (Color online) Constant energy cuts of the rotational sym- 
metrized spin dynamical structure factor in the momentum space 
in the MSW theory for 5 = 1, J1/J2 = 1.0, K/J2 = 0.8, and 
T/J2 = 1.0. The corresponding energies are respectively cu = 4J2 
in (a), LJ = IOJ2 in (b), lj = II.5J2 in (c), lj = I2J2 in (d). In all 
the panels, a broadening factor O.5J2 has been used for the conve- 
nience of calculation. 



in Ref. Isil in the antiferromagnetically ordered phase, but our 
resuhs apply to the paramagnetic phase. 



IV. DYNAMICAL STRUCTURE FACTOR 



In order to investigate the magnetic excitations, which are 
directly accessible by INS measurements, we have calculated 
the magnetic structure factor iS(q, oj). Our main interest is 
to understand the experimentally observed anisotropic feature 
of the magnetic excitations in the paramagnetic phase above 
the Neel temperature. As already discussed in Sec.|IIIl in this 
temperature regime, the most relevant factor for the in-plane 
anisotropy is the Ising order Therefore, we will concentrate 
our discussion on the magnetic structure factor in phase II of 
the 2D Ji~ J2 — K model. In this phase iS(q, w) has the same 



5(q,w) ==27r-^^' [cosh(20k+q - 2ek) - ss'] 

k s,s' = ±l 

X(5(w - s£k+q - s'ek)?^k+q?^k 7 (25) 



where refers to the summation over the magnetic Brillouin 
zone corresponding to the (tt, 0) order, which is enclosed by 

-7r/2 ^ 
and n. 



-TT 



TT. n 



+ 



kx ^ 7''/2, and 

rik. C = 1 in the MSW theory, and C 



nk + 1 

( j^ '''K" ^ ' " -L 111 LllV^ iVAvJ T* LllV^WlJ, U.11U W 3/2 111 

the Schwinger boson theory with SU(2) symmetry)^ We see 
from Eq. (l25T l that the contribution to iS(q, w) comes from 
two-boson processes. Hence, in general cases the peak of 
5(q, cj) does not follow the boson dispersion. But at low 
temperatures, the largest contribution to 5(q, uj) in the sum- 
mation over k comes from the term at k = (0,0) since the 
small gap Ai at this point results in a large boson number 7ik. 
To satisfy the energy conservation in the S function, S{q, lo) 
must be peaked at w w £q. Actually this leads to a two-peak 
structure corresponding to s' = ±1 near a; = Eq for a fixed 
q. But the separation of these two peaks is proportional to 
£(0.0) and is very small at low temperatures. In the numeri- 
cal calculations performed, the gap between the two peaks is 
healed by substituting the delta function by a Lorentzian with 
a small broadening width. As a result of this small broaden- 
ing, S{q,uj) only shows a single peak structure. Therefore, 
in this limit the peak positions of 5(q, w) follow the boson 
dispersion. 

To better discuss the anisotropic distribution of the spec- 
tral weight in momentum space, we plot the constant energy 
cuts of the calculated 5(q, uj) at a fixed temperature T < T^o 
in Fig. |4] At low energies, the peaks of 5(q, to) form a el- 
liptic ring centered at (tt, 0) (and also its symmetry related 
point (0, tt) after rotation symmetrization), as displayed in 
Figs. |4ta),(b). The elliptic feature is a consequence of the 
anisotropic correlation lengths in the Ising ordered phase, and 
the ellipticity near (tt, 0) is proportional to ^x/^y = vix/viy, 
which is not sensitive to temperature since the mean-field pa- 
rameters are only weakly temperature dependent for T < T^q 
(Fig. |2]i. The ellipticity also only weakly depends on K: for 
= 1, we find ^x/iy ^ 1-7 at X = 0, and ^x/Cy ^ 1-4 
at K/J2 = 0.8. 

With increasing energy, the ellipse centered around (tt, 0) 
expands towards the Brillouin zone boundary, as seen in 
Figs. |4ta)-(d). For sufficiently large energy, the spectral 
weight reduces greatly along the Qx direction, and the 5(q, w) 
is peaked near qy ~ ±7r/2 along the qy direction (Fig.|4jc)). 
The elliptical peak feature appears to have been split into two 
parts in the direction of its major axis. As the energy gets 
close to e(7r,7r)' the two peaks move towards (tt, ±tt), form- 
ing patterns that are centered around (tt, ±7r); cf. Figs. |4|c)- 
(d). In our theory, there are two factors that contribute to this 
anisotropic distribution of the spectral weight along the el- 
lipses. Firstly, for lj > £(,r/2.o)' along the qx axis the energy 
conservation in the S function of Eq. ( |25] | can only be satisfied 
when k ^ (0, 0). A nonzero k corresponds to a smaller nk, 
which greatly reduces S{qx,Lj)- Along the qy axis, however, 
S{qy,Ld) is not reduced because the k = (0, 0) mode can still 
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satisfy the energy conservation. Secondly, for a given k, the 
coherence factor cosh(26'k — 26'k_|_q) along the ellipse is also 
anisotropic. To see this, recall that the largest contribution to 
5(q, u) is from the k = (0, 0) term in Eq. dZST l. For simplic- 
ity, we take a single mode approximation, namely, S{q,uj) 
can be approximated by this k = (0, 0) term. Then the ellipse 
showing spectral weight peaks is determined by Eq = w, and 
the coherence factor cosh(20k—2^^k+q) oc (Sq— Aq)/eq. For 
^JqUl + <liey « r/Ai,Bq-Aq « A - Eq/ A - 2 Ji, /, 

where A = 8 J2.9x+y + 4:Jix9x- Since Jiy < for the choice 
of model parameters, it is easy to see that along the ellipse 
Eq = uj, the maximum of the coherence factor is located along 
the Qy axis but not the qx axis. Since within the single mode 
approximation, 5(q, a;) is proportional to the coherence fac- 
tor, 5(q, uj) is also anisotropic along the ellipse. Note that at 

low energies qlQ + qlQ < T/Ai), Jiyfyql < A, so 
the anisotropy is very small. This coherence-factor-induced 
anisotropy becomes sizable when the ellipse is large (for 

V. DISCUSSIONS 
A. The effects of spin size 



and found the mean-field phase diagram is similar to the one 
in Fig.[T] Approximately, r^o and Tq are increased by a factor 
of S{S + l)/2. The boson dispersion shown in Fig.jSja) also 
exhibits the similar features as in the 5=1 case. In Fig.|3b) 
we compare the ratio of the effective nearest neighbor cou- 
plings Jiy/Jix, defined in Eq. (|9]l, for 5 = 1 and 5 = 2 (at 
zero temperature). We find that with increasing 5, the mini- 
mal K/ J2 value where Jiy becomes ferromagnetic is dropped 
from K/ J2 ~ 0.53 to K/ J2 ~ 0.13. Hence, we conclude that 
the anisotropy of the effective exchange couplings induced by 
non-Heisenberg coupling K is more significant for larger spin 
size 5. 

We can further compare our MSW result at T = with the 
one in a recent MSW study, which used a mean-field treatment 
that is different from ours."*^ The two theories yield exactly the 
same results when the spin size S ~^ 00. For finite spin sizes, 
by comparing the behavior of J\yJJ\x ratio in Fig. |5lb) and 
the corresponding results in Ref. \i% we observe that the two 
theories give qualitatively similar results for the anisotropy in 
the exchange couplings: The biquadratic coupling K reduces 
the ratio of the effective ratio J\yl J\x- Quantitatively, there 
are some differences between the two approaches. In particu- 
lar, while for 5 > 1, the ratio J\yl J\x changes sign at a finite 
K value in both theories, this sign change does not appear for 
5 = lin Ref.m 



B. Generalizations of tlie J 1 — J2 — K model 
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FIG. 5. (Color online) (a): MSW dispersion along high symmetry di- 
rections in the paramagnetic Brillouin zone for the 2D Ji — J2 — K 
model at 5" = 2, J1/J2 = 1.0, T/J2 = 1.0, and various K val- 
ues, (b): The ratio of Jiy/ Jix, showing the anisotropy in effective 
exchange couplings, see Eq. l[9]l, as a function of the biquadratic cou- 
pling for 5 = 1 and 5* = 2. 

Besides the 5=1 results shown in Sec.HIIland Sec.|IV] we 
have also studied the J1 — J2 — K model with larger spin sizes. 



Several remarks on the Ji — J2 — K model studied in this 
paper. From the incipient Mott picture, when the system is 
in the vicinity of Uc/t, the spin Hamiltonian contains interac- 
tions involving more than just two spins. To see this, we start 
from a multi-orbital Hubbard model on a square lattice, and 
assume that Hund's rule coupling locks the spins in different 
orbitals to a high spin state. Then we may obtain a spin-only 
Hamiltonian by integrating out the fermion degrees of free- 
dom based on perturbation in t/U. To the t'^ /U order, we ob- 
tain the usual Ji — J2 Heisenberg interaction between nearest- 
and next-nearest- neighbor spins. The next-order terms appear 
in the order f^/U^, and include the biquadratic K term as well 
as the ring exchange interactions. Here, we have focused on 
the effects of the biquadratic interaction. The influence of the 
ring exchange interactions in the regime we are considering is 
briefly discussed in Appendix |B] 

To fully understand the antiferromagnetic (tt, 0, tt) order re- 
vealed in the experiments, the 2D Ji — J2 — X model needs to 
be extended to the 3D case by including an interlayer coupling 
Jz . A nonzero will support the antiferromagnetic order up 
to the Neel temperature T/y. In the mean-field treatment, the 
antiferromagnetic order emerges at a mean-field Neel temper- 
ature T/vo. The details of the effects of the interlayer coupling 
Jz to the magnetic phase diagram of the Ji — J2 — K model 
and the magnetic excitation spectrum is further discussed in 
Appendix ICl 

When fluctuations beyond the mean-field level are taken 
into account, the actual Neel and Ising transition tempera- 
tures, T/v and can be well below their mean-field values. 
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The mean-field temperatures T^o and Tcro then correspond to 
some crossover temperature scales, below which the fluctuat- 
ing order have significant effects. The fluctuating anisotropic 
effects we have presented will be dominant in the temperature 
regime < T < T^o^ 

C. Effect of itinerant electrons 

Within the bad-metal description of the iron pnictides, the 
quasi-localized moments are coupled to itinerant electrons 
with a spectral weight that depends on the proximity to the 
Mott transition. A convenient way to describe the effect of 
itinerant electrons on the spin dynamics is to reformulate the 
results of the local-moment-based calculations in terms of a 
non-linear sigma model, and introduce into the latter a damp- 
ing caused by the itinerant electrons; for details, we refer to 
Ref. m Well below T^yo and in the vicinity of (tt, 0), the ef- 
fects of the itinerant electrons are described in terms of the 
effective action for the staggered magnetization M: 

S[M] - T /" dq^ [r + Ar + vj^ql + vf^ql + wf 
•' I 

+-f\uji\]M^ + uM^ + 0{M'^). (26) 

Here, M = m + m is the sum of m and m', the 0(3) vectors 
respectively for the magnetizations of the two decoupled sub- 
lattices on the square lattice, and oji the Matsubara frequency. 
This action arises in a "w-expansion", which is based on a 
proximity to the Mott transition and is described in Refs.[l2l 
and [HI it has the form of the usual cr-model^. In the first 
term, Ar > is a mass shift and 7 describes the strength of 
spin damping from coupling to fermions. (See Fig.|6]l At rel- 
atively low energies, this introduces a procedure that can be 
used to describe the broadening of the spin spectral peaks in 
momentum space due to coupling to itinerant electrons^ 




FIG. 6. (Color online) Diagramm of the second-order contribution to 
the effective action in Eq. l l26t due to coupling to fermions. 

We should emphasize that this procedure is a qualitative 
treatment of the spin damping. Incorporating the full details 
of the electronic bandstructure will introduce momentum- 
dependence of the damping rate, making it possible to gen- 
erate the type of anisotropic damping that was proposed phe- 
nomenologically by Harriger et ali^. 

Comparing our results for the Ji — J2 — K model in 
Fig.lla)-(d) with those of the J1-J2 model (Fig. 4 of RefUl) 
shows that, the biquadratic term itself brings out an anisotropy 
in the spectral weight of the elliptic peaks. The spectral weight 
is larger along the major axis of the ellipse than along its mi- 
nor axis. This anisotropy goes in the same direction as that 
of the experimental data on BaFe2As2, illustrated in Fig.|4jf). 



We therefore conclude that both the ellipticity and intensity 
anisotropy of the spectral peaks in momentum space are con- 
trolled by the exchange interactions. 

We note that the Ising order parameter is also coupled to 
the itinerant electrons. Since the Ising order parameter breaks 
the symmetry, it couples to those spin singlet fermion 
bilinears that correspond to the Big representation. Conse- 
quently a nonzero Ising order parameter will induce a nonzero 
d^2_y2 nematic charge density (cos ~ cos ky)cl^ a^k,a for 
all the d-orbital electrons, where a is the orbital index. In ad- 
dition the Ising order parameter will induce a nonzero charge 

density imbalance c\.j.zC]ixz — c[^yz'^^v^ between the dxz and 
the dyz orbitals, which is also referred as the ferro-orbital or- 
der. As a result the spin fluctuations from the incoherent de- 
grees of freedom can give rise to an orbitally ordered, charge 
nematic metal, with anisotropic transport properties. In a 
model with 3D coupling (see Appendix ICb. the coupling to 
the itinerant electrons will reduce the Neel transition temper- 
ature from its mean-field value T^o to Tm < Tatq, through 
the positive Ar noted above. It will likewise decrease the 
Ising transition temperature from its mean-field value T^o to 
< Ta-o- However, the correlation lengths are still size- 
able and should be anisotropic up to To-oj^ This implies that, 
in the 3D model with three-dimensional coupling, we expect 
anisotropic magnetic excitations to exist from Tjy all the way 
up to the crossover temperature scale T^o, in the absence of a 
static Ising order 



D. Comparison witli otlier approaches 

Our studies in the Ji — J2 — K model, with or without 
the coupling to the itinerant electrons, are very different from 
purely itinerant studies with ?7/< much smaller than Uc/t. Be- 
cause the Fermi surface comprises small electron and hole 
pockets, such calculations are expected to yield very small 
spin spectral weight. Experimentally, the total spectral weight 
is known to be large, with an effective moment that is larger 
than 1 yUs/Fe in CaFe2As2 (Ref. I21I). Such a large spectral 
weight arises naturally in our approach using as the starting 
point the Ji — J2 model (with or without the K term). 

Our approach can be compared more closely with that of 
the DMFT studies of Ref. 41, in which the ratio of the effec- 
tive interaction (combined Coulomb and Hund's interactions) 
to the characteristic bandwidth is close to the Mott-transition 
value, U/t ^ Uc/t. The proximity to the Mott transition en- 
sures that a large part of the electronic spectral weight lies 
in the incoherent regime, which will naturally give rise to a 
large spin spectral weight. The consistency of the momentum- 
dependence determined by the DMFT calculations and that of 
our Ji — J2 — K calculations further suggests the compati- 
bility of the two approaches. There are some important dif- 
ferences, however. In the DMFT calculation, the anisotropy 
of the structure factor has been attributed to the geometry of 
the Fermi surface(s). The Ji — J2 — K results however tie 
the anisotropy of the spin spectral weight in momentum space 
with the Ising correlations. 

Experimentally, the Ising correlations can be very naturally 
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connected with the x~y anisotropy observed in ARPES^ and 
transport measurements in the detwinned 122 iron pnictides 
at temperatures above T^r. A recent theoretical calculation-'* 
shows how resistivity anisotropy in the tetragonal phase above 
Tjv follows from the existence of the Ising correlations dis- 
cussed here. 



VI. COMPARISON WITH EXPERIMENTS ON THE 
PARAMAGNETIC PHASES OF PARENT 122 IRON 
PNICTIDES 

Spin dynamics in the paramagnetic phase of the par- 
ent 122 iron pnictides has been recently studied via INS 
measurements For CaFe2As2, spin dynamics at low en- 

ergies (below 70 meV) has been studied by Diallo et ali^ It 
is found that the peaks of the dynamical structure factor form 
anisotropic elliptic features at low energies, similar to the re- 
sults in the antiferromagnetic phase.-' More recently, Harriger 
et al— measured the spin dynamics of BaFe2As2 up to 200 
meV. At low energies, they found the distribution of spectral 
weights in the momentum space forms similar elliptic feature 
as in the CaFe2As2 case. With increasing energy, the elliptic 
feature expands towards the Brillouin zone boundary. More- 
over, they determined the magnetic dispersion to be peaked 
(or flat-topped) near (tt, tt). Similar results have also been re- 
ported for SrFe2As2.-i2 

Our study on the Ji — J2 — K model have already pro- 
vided valuable information for understanding these experi- 
mental observations. In real materials, the various fluctua- 
tion mechanisms and the coupling to fermions/phonons will 
reduce the Neel and Ising transition temperatures. However, 
below the mean-field Ising temperature T^vi, the effective cou- 
plings between the nearest neighbors are always anisotropic. 
Hence we expect the magnetic fluctuations to be anisotropic 
for Tm{< Tc) < T < TcrQ, which corresponds to the upper 
portion of region II in Fig. [T] This anisotropy is reflected in 
the spin dynamics in the paramagnetic phase. 

To be specific, the anisotropic elliptic feature at low ener- 
gies observed in CaFe2As2 and other parent 122 compounds 
can already been understood within the Ji — J2 model.-- We 
have shown in Fig.|4]that the J1 — J2 — K model gives the simi- 
lar low-energy elliptic feature. It will be important to measure 
the spin dynamics at high energies in this material. 

Our calculated evolution of this elliptic feature as the en- 
ergy is raised in the Ji — J2 — K model can be systematically 
compared with the experimental observations in BaFe2As2 
and SrFe2As2. To see this, we fit the peak positions of cal- 
culated 5(q, oj) to the experimental magnetic excitation dis- 
persion data in BaFe2 As2, from which we can extract the best 
fitted values of the exchange couplings. Assuming 5* = 1, 
we find the fitted exchange couplings are J2 = 17 ± 4 meV, 
J1/J2 = 1.0±0.5,and A7J2 = 0.8±0.3. We find that a very 
broad range of the Ji / J 2 ratio can all fit the experimental data 
quite well. As illustrated in Fig. [3jb), any dispersion curve 
within the shaded region fits the experimental data within er- 
ror bars. But to fit the dispersion data near the local maximum 
at (tt, tt), a moderate K/ J2 ratio is necessary. For = 1, we 



find K/ J2 « 0.8 fits the data the best. On the other hand, for 
S = 2, the best fitted K/ J2 ratio is substantially reduced to 
about 0.2. 

For BaFe2As2, detailed measurements in the momentum 
space have been reported by Harriger et alM. This allows us 
to see that the agreement between our theory and the experi- 
ment is not only for the dispersion, but also for the anisotropic 
distribution of the spectral weight of iS(q, lo) in momentum 
space. 

In order to make a comparison with experimental data, we 
use Eq. ( |25l ) in the calculation of 5(q, cj) and approximate the 
delta function by the following Lorentzian broadening 



5{oj - Ae) 
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TT (CJ - Ae)2 + 72 



(27) 



Here we have assumed that the broadening mainly comes 
from the damping effect due to coupling to itinerant electrons. 
It is then reasonable to take the phenomenological broaden- 
ing factor to be the damping 7 introduced in Eq. (|26] | since in 
either the MSW or Schwinger boson theory, the damping is 
still due to the same bubble in Fig. |6] Calculating the mag- 
nitude of 7 requires a detailed microscopic theory and is be- 
yond the scope of this article, however we can use Ref. [3^ 
for reference, where it has been determined that 7/ J2 ~ 3 
for CaFe2As2. Here we assume that this ratio still holds for 
BaFe2As2 and the damping is isotropic. In Figs. |7|a)-(d) we 
replot the theoretical dynamical spin structure factor in Fig.|4] 
with this damping factor, and compare them with the experi- 
mental data in Ref. [s^ At low energies, our theory correctly 
captures the elliptical feature centered at (tt, 0) as displayed in 
Figs. 13 a), (b). Experimentally, this is seen as a filled elliptical 
spot due to damping effect, which is also shown in our theo- 
retical plot in Fig.lTJa). The evolution of the elliptical feature 
with increasing energy is also consistent with the experimen- 
tal observation: as the ellipse expands towards zone boundary, 
it gradually splits into two parts, and forms a pattern around 
(tt, tt) (see Figs. [Tjc), (d), and (f)). We reiterate that such 
anisotropic features are the properties of our Ji — J2 — K 
model either with an isotropic or without additional damping 
due to itinerant electrons. While anisotropic damping pro- 
posed in Ref. [s^ could reinforce the effect, it is not necessary 
to understand the INS experiments. In CaFe2As2 the ellipti- 
cal feature around (tt, 0) persists up to high energies, while in 
BaFe2As2, this elliptical feature splits into two parts at inter- 
mediate energyi^ These two different behaviors can both be 
understood within our Ji — J2 — K model with similar, nearly 
isotropic damping but different K values. 



VH. CONCLUSIONS 

In this paper we have investigated the finite temperature 
spin dynamics of a Ji — J2 — K antiferromagnetic Heisen- 
berg model using both MSW and SBMF theories. The spin 
dynamics obtained from these two methods are similar to each 
other 

We have found that by including a moderate biquadratic 
coupling K, the magnetic excitation spectrum of the Ji — 



10 





FIG. 7. (Color online) Evolution of 5(q, oj) in the paramagnetic phase of the Ji — J2 — K model, showing that the elliptical features near 
(tt, 0) at low energies (top panels) are split into features that are centered around {n,n), as the energy is increased towards the zone-boundary 
spin-excitation energy (bottom panels). This trend is consistent with the inelastic neutron scattering experiments, shown in the box on the 
right for two energies measured in the paramagnetic phase of BaFe2As2 (data talcen from Ref.ls^. (a)-(d): Same as Fig. |3a)-(d), but with 
damping 7 = 3 J2. (e)-(f): The INS data at T = 150 K in BaFe2As2, taken from Ref.ls^. The energy transfer is a; = 50 ± 10 meV in (e), and 
Lo = 150 ± 10 meV in (f). Here we find that the best agreement between theory and experimental data achieves when taking J2 w 13 meV in 
the model. 



J2 — K model is anisotropic below a mean-field Ising tran- 
sition temperature T^o- As in the case of the Ji — J2 model"*-, 
the peak of the low-temperature dynamical structure factor 
5(q, cj) contains elliptical features near (tt, 0) in the param- 
agnetic Brillouin zone at low excitation energies. However, 
unlike the pure Ji — J2 model, the spectral intensity also dis- 
plays anisotropy along the ellipse, with the intensity being 
higher along the major axis than that along the minor axis. 
This spectral anisotropy accounts for the observed particu- 
lar way in which the low-energy elliptical features, centered 
around (tt, 0), expand towards the zone boundary as the en- 
ergy is increased towards the zone-boundary spin-excitation 
energy. It also gives rise to a particular form of high-energy 
spectral features that are centered around (tt, tt). 

We have also compared our calculated dynamical spin 
structure factor of the Ji — J2 — K model with the recent in- 
elastic neutron-scattering measurements in the paramagnetic 
phases of the 122 iron pnictides^2i^2il£. The theoretical results 
provide a very natural understanding of the salient features of 



the experiments. 
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Appendix A: Ising transition at small Ji / J2 ratios 

We find that the nature of the mean-field Ising transition at 
Tcto depends on both J1/J2 and K/ J2 ratios. At K ~ and 
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J\l Ji ^ 0.9, we find T^o < To, and the Ising transition at 
Tjro is always second-order (Fig. [8]l. When J1/./2 ^ 0.9, T^q 
meets Tq and the Ising transition becomes first-order. This is 
an artifact of the mean-field approximation since the transition 
at To is always first-orderi^ Still for J1/J2 ^5 0.9, increas- 
ing K from zero, the transition at Tcro changes from second- 
order to first-order when K is bigger than a bicritical point 
K^. As shown in Fig.[8]for J^jJi = 0.6, A'c/ J2 ~ 0.04. At 
K >j Kc, Tcro < Tq. This suggests that the Ising transition 
near K^. is not influenced by Tq , but the order of this transi- 
tion is tuned by K. Hence the first-order transition at Tg-o is 
not an artifact of the mean-field treatment. 




0.0 0.2 0.4 



FIG. 8. (Color online) Mean-field magnetic phase diagram in the 
MSW theory for S = 1, J1/J2 = 0.6. The dashed blue and dashed 
dotted brown curves refer to the mean-field temperature scales Tno 
and To, respectively. The thicker solid red curve refers to a second- 
order Ising transition at T^o, while the thinner solid red curve refers 
to a first-order transition. In the shaded region, the effective exchange 
coupling Jiy < 0. 



Appendix B: Effects of ring exchange couplings 

Besides the quadratic and biquadratic interactions, other in- 
teractions involving more than two spins can also appear in 
the spin Hamiltonian in the vicinity of Mott transition. For 
instance, the four-spin ring exchange interaction can appear 
as a consequence of the fourth-order perturbation associated 
with the electron hopping process. We can consider the ef- 
fects of a four-spin ring exchange process on the spin dy- 
namics by adding a term A'n X^jj.a..; [(Sj • S_,)(Sfc • S;) - 
(S, • Sfc)(Sj • Si) + (Sj • S,)(Sj '• Sfc)] to the Hamiltonian, 
where TTq > 0, and the sites (i^j, k, I) are the vertices of a 
square plaquette, labeled clockwise. The four spin ring ex- 
change competes against Ji and J2 and tends to weaken the 
antiferromagnetic order coming from Ji or J2. In the linear 
spin wave description of the (tt, 0) ordered state, we obtain 
the effective exchange constants Jix = Ji + 2{K — Kij)S'^, 
Jiy = Ji- 2{K - Ka)S^, and J2 = J2 + KqS^, and a 
reduced spin gap at (tt, tt). This trend also persists in the para- 



magnetic state, and reduces the size of the Ising order param- 
eter. For consistency with the experimental results we require 
K > An. 



Appendix C: Effects of interlayer exchange coupling 
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FIG. 9. (Color online) Mean-field magnetic phase diagram in the 
MSW theory for 5 = 1, J1/J2 = 1, and an interlayer exchange 
coupling Jzl J2 ~ 0.1. The dashed blue and solid red curves refer to 
the mean-field temperature scales Tnu and T„o, respectively. In the 
shaded region, the effective exchange coupling Jiy < 0. 




FIG. 10. (Color online) The temperature evolution of the mean-field 
parameters in the MSW theory for S = 1, J1/J2 = 1, A'/J2 = 0.8, 
and with an interlayer exchange coupling Jzj Ji ~ 0.1. 
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FIG. 1 1 . (Color online) Dispersion of the J\ — J2 — K model in 
the MSW theory for various K values at S = 1, J1/J2 = 1.0, 
T / J2 = 0.1, and with an interlay er exchange coupling Jzj Ji = 0.1. 



The real materials have a 3D tetragonal structure. In the 
J\ — J2 — K model, the 3D effects can be studied by extend- 
ing the model to include a finite interlayer exchange interac- 
tion Jz Si • Si+£. In 3D the long-range antiferromagnetic 
phase survives at finite temperature up to the Neel tempera- 
ture Tn- In MSW and SBMF theories, the mean-field Neel 
temperature Tjvo is determined by the onset of spontaneous 
sublattice magnetization toq- In general, T^o ^ T„ Tq. 
The modification to our discussion in Sec. HIl comes through 
an additional interlayer antiferromanetic bond correlation pa- 
rameter gz- In the presence of , the self-consistent equations 
of Eqs. (|4li-(|6]l and Eqs. (fT8Tl-(l20li are unchanged, but the ex- 
pressions for and i?k are modified according to 



= Ak + 2J,g,cosfc, (CI) 
= -Bk + 2J,53, (C2) 



in the MSW mean-field theory, and 

= A + 2J,.g,sinfc, (C3) 
B^D = B^. (C4) 

in the SBMF theory. 

In Fig.|9]we show the phase diagram at the experimentally 
suggested ratio Jz/'h = 0.1. Similar to the 2D case, the 
mean-field phase diagram consists of an Ising and Neel or- 
dered antiferromagnetic phase (I), an Ising ordered but Neel 
disordered paramagnetic phase (II), and an Ising and Neel 
disordered paramagnetic phase (III), separated by mean-field 
temperatures Tjvo and T^q (see also Fig.fTOli. For the param- 
eters in Fig. |9] the transitions are both first-order, and both 
Tjvo and To-o increase with K. For K/J2 ^0.2, T^vo meets 
Ta-o, and there is only a single transition between phases I and 
III. The absence of phase II in this regime is an artifact of the 
mean-field theory, since Ta-o is always bounded above by the 
mean-field scale Tq- 

In connection to the real materials, we note that the struc- 
tural and magnetic transitions in the 1111 pnictides are well 
separated. But in 122 compounds, they are either very close to 
each other, or become a single first-order transition. This can 
be understood in terms of the present theory, provided Jz is 
stronger in the 122 materials. By comparing Fig.|9]and Fig.|2] 
we see that the magnetic transition is closer to the Ising tran- 
sition for a larger Jz- Recent experiments also show that the 
electron doping may cause the separation of the structural and 
magnetic transition temperatures in Ba(Fe,Co)2As2 systemi^ 
The similarity between this behavior and the K dependence 
of Tcto and Tjvo in the phase diagram of Fig. |9] suggests the 
possibility that electron doping is positively correlated with a 
reduction of the biquadratic interaction. It would then be inter- 
esting to reveal the link between them in future experimental 
and theoretical studies. 

In Fig. [TT] we show the low-temperature boson disper- 
sions of the 3D model for various K values along two high- 
symmetry directions in the fc^ = tt plane. Aside from a larger 
gap at (0, 0, tt), the dispersion is very similar to the one in 2D: 
the dispersion is highly anisotropic, and with increasing K, 
the local minimum at (tt, tt, tt) turns to a maximum. This is 
not too surprising because the in-plane anisotropy is a conse- 
quence of the 2D Ising-type fluctuations, and is not sensitive 
to the interlayer exchange coupling. 
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